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The recent uncovering of the Fermi Bubbles/haze in the Fermi 7-ray data has generated theo- 
retical work to explain such a signal of hard 7-rays in combination with the WMAP haze signal. 
Many of these theoretical models can have distinctively different implications with regards to the 
production of high energy neutrinos. We discuss the neutrino signals from different models proposed 
for the explanation of the Fermi Bubbles/haze. More explicitly, from Dark Matter annihilation in 
the galactic halo with conditions of preferential CR diffusion, from recent AGN jet activity, from 
periodic diffusive shock acceleration, from stochastic 2nd order Fermi acceleration and from long 
time-scale star formation in the galactic center in combination with strong galactic winds. We find 
that some of these models will be probed by the IceCube DeepCore detector. Moreover, with a km 3 
telescope located at the north hemisphere, we will be able to discriminate between the hadronic, 
leptonic and the DM models. Additionally using the reconstructed neutrino spectra we will probe 
annihilation of TeV scale dark matter towards the galactic center. 
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I. INTRODUCTION 

Using the first year of Fermi LAT 7-ray full sky data, 
[l| revealed the presence of a diffuse component towards 
the galactic center that extends up to 50° in latitude. 
That component has a spectrum significantly harder than 
elsewhere in the Galaxy [1-3] and a morphology elon- 
gated in latitude to longitude, that depending on the 
exact template analysis used for its extraction from the 
full sky data, is either well defined within 2 bubbles with 
distinct edges both at high and low latitudes |2[ known 
as the "Fermi Bubbles", or is slightly more diffuse (the 
"Fermi haze") with a latitude to longitude axis ratio of 
~ 2 [3] but still confined within | I \ < 20° and | b | < 50° 
and with the "edges" seen only at the higher latitudes 

a. 

That signal of hard 7-rays together with the WMAP 
haze 0, Hf may indicate the presence of a hard compo- 
nent of Cosmic Ray (CR) electrons, that through their 
up-scattering of low energy photons produce the hard 
spectrum of (Inverse Compton) 7-rays and through their 
synchrotron radiation the hard observed spectrum at mi- 
crowaves. Various authors have suggested mechanisms 
for the origin of these CR electrons. Among them, pos- 
sible scenarios include recent (1-3 Myr ago) AGN jet ac- 
tivity in the galactic center p, 0] ; TeV scale Dark Mat- 
ter (DM) annihilating to leptons Q , within conditions 
of preferential diffusion perpendicular to the plane 
stochastic 2nd order Fermi acceleration by large scale 
turbulence in magneto-sonic waves [9[, or from periodic 
injection of hot plasma causing diffusive shock accelera- 
tion (1st order Fermi acc.) in the halo fioj . 

Alternatively, CR protons associated with long time 



scale (~ Gyr) star formation in the galactic center, trans- 
ferred by strong winds into the Fermi Bubbles region 
have been suggested by [nj|. Finally, a combination of 
DM and millisecond pulsars in the galactic halo adding 
up to the signals at 7-rays and microwaves has been dis- 
cussed in |12| . 

The detection or lack of high energy neutrinos from 
km 3 neutrino telescopes could help discriminate between 
the leptonic 0, El, S S EqI and the hadronic [ll[ scenarios, 
since the leptonic scenarios would not produce any, or a 
few and up to the TeV scale neutrinos while the hadronic 
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would produce abundant neutrinos up 
lH, EH (see also discusion in [HI). 
Additionally to the discrimination among different 
models for the Fermi haze/Bubbles via searching for 
their neutrino counterpart, such searches can be used as 
an other channel of indirect DM searches. 

DM composes approximately 85% of the matter den- 
sity of the universe, yet its particle physics properties 
still remain unknown. Measurements of cosmic ra ys ( CR) 
[I3422I have generated new model building [8, 23M34j and 
have helped place new constraints on dark matter proper- 
ties [35|439[. Since many of these models and constraints 
are placed at the TeV mass scale and suggest /refer to 
enhanced (boosted) annihilation cross-sections with hard 
spectra for the SM particle annihilation products, neu- 
trino signals from the galactic center (GC) [40] from the 
Sun [41] or the Earth [42] can be considered of interest. 

IceCube in the South Pole [H, Q has already pre- 
sented some early (pre-DeepCore) results |45l-[47j and is 
expected with its DeepCore update to probe better the 
sensitive for DM searches region of ^100 GeV in neutrino 
energy. ANTARES located in the Mediterranean is also 
collecting data [48-50] and due to its location probes sig- 
nificantly better the GC than the IceCube, but because 
of its small size has still low number of statistics. Fi- 
nally a future km 3 telescope located in the Mediterranean 
as KM3NeT [111 (EH, [Bl will combine the virtues of 
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IceCube and ANTARES providing a good neutrino tele- 
scope for searches of DM annihilation towards the GC. 

In this paper, in section [II] we will discuss the neu- 
trino signals of various DM models that could explain the 
Fermi haze and WMAP haze signals as presented in [3|, 
for both IceCube and a future km 3 telescope located in 
the Mediterranean, presenting also search strategies for 
those signals. We use simulated performance information 
published for the KM3NeT [EH, E2). We will extend in 
section IIIII the discussion of searching for signals in neu- 
trinos towards the GC from more generic DM models. 
In section HVl we will describe the neutrino predictions of 
the hadronic model of [11] for the Fermi Bubbles that if 
true should soon be seen. In section [V] we discuss why 
we don't expect any significant signal in neutrinos from 
the non-DM leptonic models of [6, 7, 9, 10] presented for 
the Bubbles and conclude in section [VTl 

In this work we will discuss only the upward going 
v^s flux expected to be measured at the ongoing and 
future telescopes, as contained events numbers are ex- 
pected to be significantly suppressed by comparison; see 
for instance discussion in [53j. 



II. EXCITING DARK MATTER 
ANNIHILATION, CONNECTING TO THE 
FERMI HAZE AND THE WMAP HAZE 

In the context of DM annihilating to leptons, Q in- 
voked preferential (anisotropic) diffusion of CR electrons 
perpendicular to the galactic disk due to ordered mag- 
netic fields in the same direction. Such conditions could 
explain the Fermi haze spectrum and morphology in 
combination with the WMAP haze spectrum and angu- 
lar profile (see Figs. 7 and 8 of |3[). 

In that case, eXciting Dark Matter (XDM) was 
considered where DM particles with mass m x = 1.2 TeV 
annihilate into a pair of scalar bosons <fi that then decay 
to a pair of e ± due to kinematic suppression for the case 
of < 2m M . In such a case the energy released by 
the DM annihilation goes to e ± that have a hard enough 
spectrum to explain the haze signals [3|, [8[ . 

Also in an XDM scenario as has been shown in [23] , the 
annihilation cross-section can be enhanced up to O(10 3 ) 
to motivate the necessary boost (BF) on the annihilation 
rates that have been invoked to explain the CR positron 
fraction and e~ + e + flux excesses Thus the 

magnitude of the annihilation rate needed to explain the 
Fermi and WMAP haze signals (BF~ 30), is naturally 
explained within the context of an XDM annihilating to 
final state SM leptons. Finally, since from Cold Dark 
Matter (CDM) cosmological simulations [57l459| the DM 
halo profiles are typically triaxial, a prolate DM profile 
has been used in j3j. 

For the case when the SM leptons from the <p decay 
are e^ (m^ < 2m M ), no neutrinos are produced thus the 
presence of a possible "Neutrino haze" is excluded. Al- 
ternatively, for > 2m M the decay to p ± is allowed 



and a neutrino haze can exist. To maximize the possible 
neutrino haze signal and also to explain the Fermi and 
WMAP signals we will consider the case where m x = 2.5 
TeV particles an nihil ate to a pair of 0s that decay with 
a BR=1 to /i ± fl20) . The mass of 2.5 TeV is chosen 
to produce -after the muons decay- e ± which during 
propagation will give similar synchrotron and IC signals. 
The necessary enhancement in the annihilation rate is 
BF~150 to produce the same total injected energy in 
high energy e ± (see for instance Fig. 6 & 7 of p55|). 

Since neutrinos do not diffuse or loose energy, the neu- 
trino signal from the DM halo will be identical to the 
annihilation rate profile ~ J (av) p 2 DM dldVt] with / being 
the line of sight, dVt the angle of observation, p^M the 
DM density and (av) the velocity averaged annihilation 
cross-section (more accurately written as (a | v |)). Since 
in the Sommerfeld enhancement case, the (av) depends 
on the velocity dispersion [23|, [60l - [62| , it also may have 
a profile within the main halo. That is ignoring effects 
of substructure that may make the position dependence 
of the averaged (after integration) (av) over the galaxy 
even more evident 0). 

In Fig. [1] we show the case of 3 x 10 4 neutrino simulated 
events for a prolate DM Einasto profile described by: 
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with p = 0.4 GeV cm -3 the local DM density [H, |gg|, 
a = 0.17, Z c /R c = 2 and Z c = 27 kpc (giving a total 
amount of DM within inner ~ 100 kpc as is for a spheri- 
cally symmetric case of Z c = R c = 25 kpc) [66[. 

The cases of a homogeneous (av) Fig. [T] (left) and of 
(av) oc r 1 / 4 in Fig. [T] (right), (with r the galactocentric 
distance) are shown. The latter dependence can be the 
case for the Sommerfeld models [8j , as cosmological sim- 
ulations with DM and baryons have suggested a profile of 
increasing velocity dispersion towards the GC (compared 
to the local values) and up to the inner 1 kpc |68l-[71j. 
The profile with the ex r 1 / 4 cross section is slightly more 
diffuse and less elongated. Since the exact DM profile in 
the inner 5° of the galactic center is very uncertain we 
will derive our conclusions ignoring that part of the DM 
halo. Additionally, TeV neutrino sources either point or 
diffuse (from inelastic collisions of CR nuclei with dense 
ISM gas) are concentrated along the galactic disk. The 
thin galactic disk where continuous star formation takes 
place, has a characteristic (for exponentially decreasing) 
scale high of 0.3 This scale hight is indi cative 

for the majority of the TeV neutrino sources [121] . The 
HI ISM gas has also a scale hight that towards the GC is 
~ 0.15 kpc and increasing up to 0.3 kpc at the solar ring 
[lIllI; while the H2 gas a scale hight of 0.1 kpc |zSlzi|. 
These scale hight s are indicative for the galactic diffuse 
TeV neutrino flux |l22| (see also work of |77||). Thus we 
will avoid the entire inner 5° in | b |. 

The neutrino flux at Earth due to DM annihilation 
from an angle dVt ignoring oscillation (<fi®i) is described 
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FIG. 1: 3 x 10 4 simulated v events, from XDM to fi^ scenario. Left: Prolate Einasto profile with homogeneous enhancement, 
right: including a oc r 1 ^ 4 in the annihilation cross-section. The latter case is less prolate in its morphology. We present the 
neutrino maps in Mollweide projection using HEALPix [671 ]. 
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where we have left (<jv}(£, 0) to depend on the position in 
the Galaxy for the most generic case. A boost factor is 
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absorbed in either p 2 DM or/and the (crv). The is the 
neutrino spectrum of the species V{ . The multiplicity M % 
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originating from the 2.5 TeV XDM to 
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In this work we are going to discuss only the up- 
ward going z/ M s flux. The A e ff of v e upward for both 
the IceCube DeepCore (not optimally placed for the GC 
searches) [44] and the KM3NeT jzjf, is smaller for v e s 
by at least a factor of 2 at all energies of interest. For 
simplicity we are going to ignore their contribution. 

The observed flux at Earth after oscillations is given 
by [45 F 



of vis per annihilation event is absorbed in 70 ff giving: muons case, the injection spectra of v^, v e , v e (there 

are no v T s) are practically identical to those of the in- 
jected e ± given in appendix A of [81]. Per annihilation 
event there are 2 neutrinos (and 2 antineutrinos) for each 
flavor. 

Having excluded the | b |< 5° region, the basic re- 
maining background is that of the atmospheric upward 
neutrinos. The atmospheric background flux is isotropic 
after averaging for the many different directions of the 
neutrino telescopes axis within long timescales. For the 
atmospheric v ^ and spectra and fluxes above 10 GeV 
and up to 10 TeV we used the tables of Appendix B 
of [HI, extrapolating to higher energies with a spectral 
power law of 3.7 for the differential spectrum [ioj . 

In Fig. [2] we give in galactic coordinates 10 yr mock 
maps for + upward events with energy between 
360 and 2160 GeV. The energy range has been chosen to 
optimize the detection of a DM signal for the specific 2.5 
TeV XDM case. 

The atmospheric background + V ^ events that are 
shown in the top left of Fig. [2] are 152784; while the DM 
events for the prolate Einasto are 424 and 332 for the 
case of (crv} oc r 1 / 4 for the entire sky. Those numbers 
of DM events are smaller than the 1 a deviation (for the 
entire sky). Since the morphology of the DM signal is 
much different than that of the atmospheric background, 
one can expect to see a signal increase of events towards 
the GC (see the bottom row of Fig. [2|) where we show the 
inner 60° x 60°. Including the | b \ < 5° mask for the TeV 
sources (point and diffuse) for the galactic center and disk 
most of this dim DM contribution is hidden. In fact since 
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where <\P V , is the flux at injection of flavor species v\. 

For specific experiments one has to include the strong 
dependence of the telescopes effective area with angle 
and energy. Within an angle dVt and an energy range E - 
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FIG. 2: IceCube DeepCore z^x+F^ events after 10 years of data collected with "online filter" [44( with energy between 360-2160 
GeV. For the calculations we use the angular dependence for A e // of [43j]. Top left: Just atmospheric background, 152784 v^+V^ 
events. Top right: atmospheric background (152784 v^+V^ events) and DM annihilation contribution (424 v^+V^ events) from 
prolate Einasto profile with homogeneous annihilation cross-section enhancement. Middle: atmospheric background (152784 
+ events) and DM annihilation contribution (332 + F M events) from prolate Einasto profile with oc r 1//4 annihilation 
cross-section enhancement. Numbers refer to the entire sky. We mask out the | b |< 5° to account for neutrinos from point 
sources concentrated on the disk and from galactic diffuse concentrated also on the disk (see text for more details). Bottom 
left: inner 60° x 60° for the atmospheric and the XDM prolate profile case. Bottom right: same region as in bottom left for the 
case of atmospheric background and XDM prolate Einasto with radial dependence. The neutrinos from the GC are minimally 
enhanced. 



for IceCube the neutrino events are minimally enhanced 
in the GC and TeV neutrino point sources and the diffuse 
galactic neutrinos are also expected to contribute in that 
region of the sky, a claim for DM can not be made. Thus 
for the IceCube DeepCore the sensitivity towards the GC 
is still too low and the angular resolution (not accounted 
for here) too large (^ 5°) to provide a robust signal for 



that DM model. 

For a km 3 telescope in the North Hemisphere the situa- 
tion can be very different. In Fig.[3]we give the expected 3 
yr mock maps of reconstruct ed ev ents using the HOURS 
simulation [Fl| for KM3NeT |l25( . In that case by com- 
paring Fig. [3] top left where we show only the atmo- 
spheric background simulated events with Fig.[3]top right 
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FIG. 3: KM3NeT + events after 3 years of data reconstructed by the HOURS [5l|] package with energy in the range 
of 360-2160 GeV. We use for KM3NeT the same angular dependence of A eff as for ANTARES due to their expected 
similar geographic latitude location. Top left: Just atmospheric background, 212446 + events. Top right: atmospheric 
background (212446 + F M events) and DM annihilation contribution (2412 Vy, + F M events) from prolate Einasto profile 
with homogeneous annihilation cross-section enhancement. Middle: atmospheric background (212446 +F M events) and DM 
annihilation contribution (1579 v^+Vn events) from prolate Einasto profile with oc r 1 ^ 4 annihilation cross-section enhancement. 
As in Fig. [2] events numbers refer to the entire sky and we use the same mask of | b |< 5°. Bottom left and bottom right: as the 
bottom plots of Fig. [2] for the case of atmospheric background and XDM prolate Einasto without (left) and with (right) radial 
dependence. Due to its high sensitivity towards the GC and its good angular resolution KM3NeT would observe a clear signal 
from these models. 



(for combined DM Einasto prolate and the atmospheric 
background), and middle (for the prolate Einasto with 
(av) oc r 1 / 4 plus atmospheric), one can see a clear indica- 
tion of a signal from the 2.5 TeV XDM to muons model 
used to explain the combination of Fermi and WMAP 
haze signals. Concentrating in the inner 60° x 60° win- 
dow for the case of simple Einasto profile the signal will 



be very clear even in the total counts. Known TeV 7-ray 
sources such as the H.E.S.S. J1745-290 are also 

expectred to contribute in that region. Yet the angular 
resolution of KM3NeT at ~ ITeV is expected to be below 
1° [5l|, thus given the already precise known location of 
these sources (from 7-rays) it will be easy to account for 
their contribution, as has been done with 7-rays H.E.S.S. 



6 



5xl0 5 



% 
O 




5xl0 5 



5 
> 

o 



500 700 1000 

E v ^(GeV) 



1500 2000 3000 




500 700 1000 1500 2000 
E v ^(GeV) 



3000 



FIG. 4: KM3NeT +F M reconstructed spectra after 3 years within 5° <| b |< 15° and | / |< 5°. Blue solid line: atmospheric 
background flux, red bashed line: DM only flux, red solid line: atmospheric+DM flux. Left: Einasto prolate profile with 
homogeneous annihilation cross-section enhancement. Right: Einasto prolate profile with oc r 1//4 with annihilation cross-section 
enhancement. 



data [86| and with Fermi data [87j . 

At 10 yrs of data with KM3NeT the excess towards the 
GC will be very evident. About 8 x 10 3 DM events for 
the case of an Einasto prolate profile and 5 x 10 3 for the 
Einasto prolate profile with (av) oc r 1 / 4 will be detected. 

In Fig. [H we show the + reconstructed spectra 
from atmospheric and from XDM to muons after 3 yrs of 
observations in KM3NeT, from a window of 5° <| b |< 
15° and | / |< 5°, which was chosen to be optimal for a 
search of a signal from DM annihilation. 

The uncertainty in the atmospheric neutrino flux due 
to uncertainties in the p-p collisions production cross- 
sections through the decay of 7r, K and a mesons is ex- 
pected to be up to ~ 30% at TeV energies [82|. Yet, 
this uncertainty can not explain a change in the power 
law that would appear direction dependent in the galac- 
tic sky. For the case in Fig. [4] right ((cry) oc r 1 / 4 ) the 
DM signal is too small to be detected. For the case of 
the Einasto prolate of Fig. [H left, the break at ~ 2 TeV is 
not going to be very strong, even in that window. Yet, a 
gradual hardening of the total events spectra decreasing 
from high | b | towards the galactic disk (excluding the 
| b |< 5° region), will be an indication of a signal from 
DM annihilation in the main DM halo. We expect that 
the power law of total reconstructed , events with 
energy between 300 GeV and 1.5 TeV will become harder 
by 0.3 from the window of | b |> 50°, | I |< 5° to that of 
5° <| b |< 15°, | I |< 5°. 



III. DARK MATTER ANNIHILATION 
TOWARDS THE GC, OTHER CHANNELS 

Apart from the question on explaining the Fermi and 
WMAP haze signals via DM annihilation, the general 
connection of DM searches in neutrinos has received some 
attention in the recent yea rs either from annihilation to- 



The Fermi and HESS e~ + e+ CR flux [HHUI and 
PAMELA positron [16j excesses apart from the XDM 
models discussed in [II] can be explained by a variety of 
other phenomenological channels [53, l90l - [93j . 

For standard phenomenological models XX — > 
W+W~ , XX — > bb strong limits have been placed us- 
ing SUPER-K and IceCude observations towards the Sun 
(see for example some recent works of [4l|, [88], l89|). 

Among the many models we choose to show results 
for the phenomenological models of DM annihilating 
directly to muons (xx — > or directly to W 

bosons (xx — > W+W~), where electroweak corrections 
-especially important for the former channel- have been 
included [94], [95|. These two channels are distinct be- 
tween each other since for the one annihilating to muons 
high energy neutrinos come from the decay of the boosted 
muons, while for the XX — > W+W~ neutrinos are pro- 
duced with a softer overall spectrum but at significantly 
higher multiplicity. 

In Fig. [5] we show for the window of 5° <| b |< 15°, 
| I |< 5° the expected reconstructed upward fluxes in 



IceCube DeepCore of z/„ 



for XX 



1.5 TeV (left) and for X X — ► W+W~ 



with 
with ra v = 



wards the GC [40|, l4d l53l o r from the DM annihilation 



captured in the Sun [4l|, |45|, |88|, \8£ 



2.0 TeV (right). 

The masses and annihilation cross sections are chosen 
to fit the Fermi, HESS and PAMELA leptonic excesses 
given in 1 55]. We show results for 2 cases of DM pro- 
files, the Einasto prolate profile of eq. [T] (red lines) and 
the spherical Einasto profile (green lines). As can be 
seen for the more optimistic (to give a clear DM signal) 
prolate profile, a signal can be seen in IceCube Deep- 
Core within 3 yrs of data for both channels. For the less 
optimistic spherical Einasto profile only in the case of 
XX — > W+W a weak break at ~ 2 TeV may be seen. 

Yet, in KM3NeT -that is near to optimal for searching 
for such a signal- after 3 yr of data the reconstructed 
spectra in the same window will be measured with much 
greater statistics as is shown in Fig. [6l Apart from the 
case of XX — > with a spherical Einasto profile, 
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FIG. 5: IceCube DeepCore +F M reconstructed spectra after 3 years in the window of 5° <| b |< 15° and | / |< 5°. Blue solid 
line: atmospheric background flux, red dashed and green dotted lines: DM only flux, red solid and green dashed-dotted lines: 
atmospheric+DM flux. Left: XX — > Einasto prolate profile (red dashed/solid lines) , Einasto spherical profile (green 

dotted/ dashed- dotted lines) with homogeneous annihilation cross-section enhancement. Right: XX — ^ Einasto prolate 

(red dashed/solid lines) , Einasto spherical (green dotted/ dashed- dotted lines) with homogeneous annihilation cross-section 
enhancement. 




FIG. 6: KM3NeT + reconstructed spectra after 3 years in the window of 5° <| b |< 15° and | / |< 5°. Blue solid 
line: atmospheric background flux, red bashed and green dotted lines: DM only flux, red solid and green dashed-dotted lines: 
atmospheric+DM flux. Left: XX — > Einasto prolate profile (red dashed/solid lines) , Einasto spherical profile (green 

dotted/ dashed- dotted lines) with homogeneous annihilation cross-section enhancement. Right: XX — ^ Einasto prolate 

(red dashed/solid lines) , Einasto spherical (green dotted/ dashed- dotted lines) with homogeneous annihilation cross-section 
enhancement. 



all other cases give a clear break at the respective mass 
of the DM particle. For the case of XX — > in 
the spherical Einasto an indication of a signal will be the 
smooth hardening of the spectrum as one moves (using 
the same longitude range) from the high | b | towards 
the disk as discussed in section HU In the specific case at 
| b |> 40° the power law is going to be ~ 4.0 while at the 
5° <| b |< 15° it is going to be ~ 3.2 between energies 
of 300 GeV and 1.0 TeV. For a more cored profile the 
difference in the power law between the | b \ > 40° region 
and the 5° <| b |< 15° region of the sky, is going to be 
only by 0.2. 

An alternative indication of how much better a km 3 
telescope in the North Hemisphere is going to be com- 
pared to IceCube DeepCore for that type of searches, 
is given in Fig. [7| In Fig. [71 we show the reconstructed 
^/i + within 500 GeV and 1.5 TeV for the m x = 1.5 



TeV, XX — ^ case. For the IceCube DeepCore 

(left) and the KM3NeT (right) using the reconstructed 
HOURS simulation [5l|. For the entire DM halo (ignor- 
ing substructure) we show 368 + events (3 yr) in 
IceCube DeepCore and 6482 for KM3NeT (3yr), with the 
equivalent background events to be 2.2 x 10 4 and 1.4 x 10 5 . 

We note that current limits on the muon channel 
(strongest of which come from dwarf spheroidal galax- 
ies (dSph)[96|, [§3) can not exclude a cross section of 
(av) ~ 9.5 x 10 -24 cm 3 s _1 used here for the 1.5 TeV 
mass. For the case of XX — > W+W~ the cross section 
was taken to be c± 7 x 10 -23 cm 3 s _1 . This cross section 
is a factor of ^ 5 times higher that conservative limits 
coming from dSph or antiprotons jsfj I!?], HH. F° r the 
most optimistic cross sections that are still allowed by 
ps and 7-rays from dSphs, only for a very optimized DM 
halo profile (such as a prolate or in general cuspy profile), 
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FIG. 7: Vn+Vn events with energy between 500 GeV and 1.5 TeV from DM annihilation of M x = 1.5 TeV XX — > M M - Einasto 
prolate profile. Top Left: With IceCube DeepCore in 3 yr, with "Online Filter" atmospheric background (22477 v^+V^ events) 
and contribution from DM (368 + F M events). Top Right: With KM3NeT in 3yr using HOURS reconstruction technique, 
atmospheric background (138560 Vn+Vp events) and contribution from DM (6482 v^ + v^ events). As in Fig.[2]events numbers 
refer to the entire sky and we use the same mask of | b |< 5°. Bottom left and bottom right: Zooming in the 60° x 60° window 
for the IceCube (left) and KM3NeT (right) maps of the top row. Even in IceCube some excess of events is expected to be seen 
towards the GC. With KM3NeT sensitivity and angular resolution a clear signal from that model will be observed or strong 
constraints will be placed. 



can there be a signal that would be detected. Thus, the 
neutrino searches for hadronic channels are deemed not 
optimal given how strong constraints can be drawn from 
the local p flux and from 7-ray searches for these chan- 
nels. For leptonic to mainly muons channels though, p 
and 7-s provide weak constraints and neutrinos can pro- 
vide a useful alternative search channel. 



IV. HADRONIC SCENARIO: NEUTRINOS 
FROM THE FERMI BUBBLES 

As discussed in [ll| a possible explanation of the Fermi 
Bubbles signal [H is that of copious and long time-scale 
star formation in the galactic center giving CR protons 
with energies up to the PeV scale. These CR protons 
are transferred up to a distance of 10 kpc away from the 
galactic disk due to strong winds [11] . In that scenario 
the 7-rays composing the Fermi Bubbles will come from 
the decay chains of boosted mesons produced in p-p col- 
lisions. The same process will produce neutrinos with en- 



ergies up to the cut-off energy of these hard CR protons. 
Similarly these processes take place in the atmosphere 
producing the equivalent background. 

The CR protons entering the atmosphere have a signifi- 
cantly softer spectrum ( ~ 2. 6 7 ab ove 300 GeV) measured 
most recently from [99|, llOOj |l26| , compared to those re- 
sponsible for the Bubbles that have a spectrum described 
by 0: 

^dE = N E; 21 exp [-E P /E P0 ] , (7) 

with E Po the cut-off energy ~ PeV. Thus we can expect 
to see the neutri no c omponent from the Bubbles at high 
energies [13, EH, [l0lj . 

An other difference between the atmospheric neutrino 
background and the possible neutrino signal from the 
Fermi Bubbles, is that the CR protons entering the at- 
mosphere due to column densities of matter ~ kg cm -2 , 
produce extensive showers that can reach for the most 
energetic protons up to 10 10 particles at peak number 
|l02l - ll04j while for the CR protons at the Bubbles region 
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30000 events 




FIG. 8: Fermi Bubbles, 3 x 10 4 simulated events. 



one expects much lower column densities [127j . Thus the 
possible Fermi Bubbles neutrinos can not come from ex- 
tensive showers, where products (protons mainly) from 
the hadronization of the initial p-p collision would then 
hit on new target protons. Rather, the neutrinos will 
come from the decay chains of the hadronization prod- 
ucts related to a single hard p-p inelastic collision. This is 
an additional reason for the neutrinos from the Bubbles 
to have a harder spectrum than the atmospheric ones. 

For the p-p inelastic processes the neutrinos are mainly 
produced from charged pion decays. Fo r th e neutrino 
spectra we follow t he p arametrization of |l05( , that was 
based on SIBYLL |l06] simulations of p-p collisions and 
is optimal at energies above 100 GeV that we care for. 

The neutrinos coming from the Bubbles will have the 
same morphology as the 7-rays, which is relatively flat in 
longitude and latitude with clear edges [2] . While it may 
not be trivial for the CR protons and the ISM target 
material transferred with them by the galactic winds, 
to cause such a flat morphology in I and 6, at 7-rays; 
following the assumptions of for the 7-rays we will 
take the morphology of the neutrinos shown in Fig. [8] to 
be flat within the Bubbles region, with clear edges as in 
the Fermi Bubbles signal of [2[. 

As can be seen by comparing the morphology of Fig. [8] 
to those of Fig. [TJ that are for neutrinos from DM sce- 
narios which fit the Fermi and WMAP haze, the two 
types of morphologies are distinctively different. More- 
over the neutrino spectra and fluxes as we will show differ 
dramatically. In both Figs. [T] and [8] we show the same 
number of neutrino events without specifying the energy 
range or period of observation and we are not taking into 
account that any actual neutrino telescope has (or will 
have) a strong angular dependence of its sensitivity. This 
is done to "spotlight" the different morphologies of the 
possible neutrino signals. Even after taking into account 
the specific properties of neutrino telescopes that we show 
results for, the different morphologies lead to searching 
for these signals in different parts of the sky. 

The total energy stored in the CR protons in the Bub- 
bles is estimated to be ~ 10 56 erg due to an estimated 
averaged 10 39 erg/s of injected power to hard CR pro- 



tons transferred from the GC via galactic winds in the 
Fermi Bubbles regions. This process is estimated to have 
been ongoing for a timescale of multi Gyrs |ll|. These 
assumptions can result in a quasi-steady state injected 
energy from protons to 7, e ± and v of Q p ~ 3.6 x 10 38 
erg/s from 10 GeV to 1 PeV [11]. Of that power from 
approximate equal partition of energy to 7r°, 7r + and 7r~, 
1/3 goes to 7s giving the better estimated power in the 
Fermi Bubbles of c± 2 x 10 37 erg/s for 7-ray with energy 
1-100 GeV [2j], or after extrapolation of the 7-ray spec- 
trum, ^ 1.2 x 10 38 erg/s for energies between 10 GeV 
and 1 TeV. Also from the energy equipartition to the 
product pions we can estimate that the power to neutri- 
nos is ~ 2 x 10 38 erg/s inside the Bubbles. We take the 
power to neutrinos to be 1 x 10 38 erg/s, to account for 
an overestimation by a factor of 2 on the 7-ray luminos- 
ity of the bubbles by QjJ. Assuming isotropic emission 
of us inside the bubbles and a mean distance squared 
D 2 ~ R 2 sun — 8.5 2 kpc 2 from us, we can estimate their 
flux. Their observed morphology as shown in Fig. [8] is 
isotropic within the Bubbles. 

In Fig. [9] (left) we show the expected + z7 M events 
mock map, with energy between 100 TeV to 1 PeV from 
IceCube DeepCore after 10 yrs of data collection. We 
show 74 events from the Bubbles and 63 from the atmo- 
spheric background. In Fig. [9] (right) we also show what 
KM3NeT experiment with 3 yrs of reconstructed events 
would observe. 

For the Bubbles explanation of [ll[, masking out or 
not the disk, can not affect the arguments on detecting 
the signal since most of the signal neutrino events are 
significantly above/bellow the disk. 

Extragalactic point sources also can lead to an addi- 
tional isotropic neutrino component at very high ener- 
gies where the atmospheric background gets suppressed. 
Again the same process (p-p collisions) producing high 
energy (~ PeV) neutrinos will also produce 7-rays with 
a spectrum that extends to lower energies. Due to the 
high sensitivity of the Fermi LAT instrument at the 0.1- 
100 GeV range, the more likely extragalactic sources to 
produce neutrinos ha ve alr eady been detected as point 
sources in 7-rays as in |l07| . With KM3NeT angular res- 
olution of > 0.1° above PeV energies it will be very easy 
to associate even single neutrino events to known 7-ray 
point sources where p-p collisions is expected to be the 
dominant mechanism for their production (as for instance 
in star- forming galaxies). Alternatively on can mask out 
known 7-ray point sources. Thus extragalactic point 
sources contribution to the neutrino background can be 
accounted for. Finally cosmogenic neutrinos from Ultra 
High energy CR (UHECR) protons interacting with the 
CMB are already being constrained from t he equiva lent 
7-ray spectrum at the Fermi LAT energies jl08l - llll| and 
are expected to be a significant component only at ener- 
gies above the O(l) PeV range. Thus cosmogenic neutri- 
nos do not contribute in the maps of Fig. [9l that show 
energies of neutrinos with 100 TeV < E„ v < 1 PeV 
fl28l. 
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FIG. 9: Vjj, + events, with energy between 100 TeV and 1.0 PeV. Right: With IceCube DeepCore in 10 yr, with "Online 
Filter" atmospheric background (63 +F M events) and contribution from the Fermi Bubbles (74 +F M events). Left: With 
KM3NeT in 3yr using HOURS reconstruction technique, atmospheric background (465 + 77^ events) and contribution from 
Fermi Bubbles (1795 + 77^ events). 
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FIG. 10: + reconstructed spectra from Fermi Bubbles. Blue solid line: atmospheric background flux, red bashed lines: 
flux from Bubbles, red solid lines: atmospheric+Bubbles flux. Right: IceCube DeepCore with "Online filter" after 3 years of 
data in the window of 10° < b < 50° and 0° < I < 20°. Left: KM3NeT with "HOURS" reconstruction after 3 years of data in 
the window of -50° < b < -10° and -20° < I < 0°. 



The optimal region for IceCube DeepCore to search for 
a signal of Bubbles at neutrinos would be the left part 
of the north bubble 10° < b < 50°, 0° < I < 20°, for 
which in Fig. \W\ (left) we show the reconstructed fluxes 
from atmospheric and from Bubbles after 3 yrs. Just 
based on that, one should expect soon a detection of the 
Bubbles with IceCube DeepCore if the scenario of [11[ is 
correct, or alternatively setting constraints on the hard 
CR proton component inside the Bubbles. 

With KM3NeT after 3 yrs as we show in Fig. [9] (right) 
the optimal region is the part of the south bubble with 
-50° < b < -10°, -20° < I < 0°. A clear observa- 
tion of the morphology (i.e. the south and right edge) 
would be expected. In Fig. \W\ (right) we also show the 
reconstructed fluxes at which point, if observed at that 
level would provide a very good measurement also of the 
injected energy to neutrinos (providing an alternative es- 
timate of the CR protons energy). Alternatively lack of 
detection at KM3NeT should exclude the model of [ll| 
as a possibility for the Fermi Bubbles. 



V. AGN AND ALTERNATE EXPLANATIONS 
FOR THE FERMI BUBBLES 



An alternative explanation for the Fermi Bubbles/haze 
to those of sections HT1 and HVl is strong AGN jet activity 
in the Galaxy as 0, Q • The AGN case is supported by 
very recent evidence of 7-ray jets extending out from the 
galactic center, together with a 15° width and 40° — 45° 
length cocoon at the south galactic hemisphere |l!2( |l!3| . 
For the AGN scenario, the 7-ray signal is mainly from 
CR electrons with energies up to at least ~ TeV that 
up-scatter the local radiation field (mainly CMB at high 
latitudes). CR protons can not contribute much of the 
observed 7-ray signal since the ISM gas targets have a 
very low density at high distances (up to 10 kpc) above 
the disk. Given the CR energy density profiles from the 
MHD simulations of [6], it will be very difficult for CR 
protons to explain the relatively flat in latitude morphol- 
ogy of the Bubbles which suggests a limb brightening [9[ . 
Thus the scenario of [6[ may give only few neutrino events 
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at high latitudes. 

In the 1st [10| and 2nd [9| order Fermi acceleration sce- 
narios, also some protons would be accelerated at high 
energies with power-law spectral indices (for the differ- 
ential spectrum) of E~ 2 and E~ x [9j respectively. Since 
those protons would not loose their energy fast (com- 
pared to the CR electrons), the CR proton spectra would 
be homogeneous within the bubbles. Yet, in both scenar- 
ios the explanation of the 7-ray signal is entirely from the 
leptonic components. The protons (at least in those ba- 
sic scenarios) are not expected to contribute much to the 
7-ray signal; especially since in order for protons to con- 
tribute significantly to the Bubbles signal much greater 
amounts of total energy to accel erated CRs are needed 
than in the models of [3, [Tol |l29| . As an example, in the 
mechanism presented in [9[, the needed total energy in 
CR electrons above 100 MeV is ~ 10 51 erg, while for the 
hadronic model of [ll| the total energy in CR protons 
~ 10 56 erg. Thus we do not expect a significant number 
of neutrinos from the leptonic models of [H, [To| either. 



VI. CONCLUSIONS 

The recent uncovering of the Fermi Bubbles/haze plQ 
in the Fermi 7-ray data has generated theoretical work 
in explaining such a signal in combination (or not) with 
the WMAP haze signal of @, 

We have shown that for the DM explanation of the 
combined Fermi haze and WMAP haze as in [3] under 
the optimal for neutrinos annihilation channel to muons 
(XDM [Hi! annihilation to /i^) and for a prolate DM halo 
we can observe a counterpart signal with a km 3 telescope 
located in the north hemisphere at ~3 yrs of data collec- 
tion (see discussion in section HJ and Figs. [3]- [4]). IceCube 
DeepCore and ANTARES will not observe any signal for 
such a model; while for the case of XDM annihilation to 

< 2m M no neutri- 



through a very light scalar boson 



nos are produced. 

For other channels/models of DM annihilation that 
produce more neutrinos either from larger suggested 
boost factors as is XX — > for the explanation of 

the Fermi e + + e~ signal (see for example [39j), or due 
to large hadronic branching ratios as is XX — > W+W~ , 
the neutrino events are enhanced significantly. Some sig- 
nal from DM is expected even after excluding the disk 
region where non atmospheric backgrounds are concen- 
trated (see Fig. [7j). Yet neutrinos can not provide the 
strongest limits for the hadronic channels of annihilating 
DM. For leptonic to muons case the limits can be useful 
though, since they can be more robust than the pretty 
weak limits from 7s and p [IH-IH, Hzl, [H|. 

For the Fermi Bubbles explanation of [ll|, a signif- 
icant number of > 10 TeV neutrinos is estimated and 
even with IceCube DeepCore we should expect detection 
or limits (see Figs. l9l and [T0|) . Furthermore, a km 3 in the 
north hemisphere telescope will either exclude the model 
of [ll| or confirm the morphology of the Bubbles at neu- 



trinos and measure the injected to neutrinos power from 
inelastic pp collisions inside the Bubbles as we show in 
Figs. M and M 

Leptonic scenarios for the Fermi Bubbles such as those 
of [6|, l9|, ll0| would also predict some CR protons, but are 
not expected to give any significant neutrino signal, since 
the main source for the 7-rays is IC scattering from CR 
electrons. 

With the current IceCube DeepCore telescope a first 
probe to some of the models of the Fermi Bubbles/haze 
can be achieved, while with a km 3 telescope located in the 
north hemisphere, discrimination between the hadronic 
[HI, the leptonic @, H, [To| and the DM [3] cases will be 
attained. 

As this paper was being written, a new analysis of 
Fermi 7-rays has suggested the evidence of 7-rays jets 
in the Milky Way extending out to ~ 10 kpc from the 
galactic center [112| . Additionally a cocoon structure has 
be en re vealed in the south hemisphere. While according 
to |ll2j the total luminosity of the north and the south 
jet-like features is (1.8 ±0.35) x 10 35 erg/s at 1-100 GeV 
i.e. 2 orders of magnitude less in luminosity than the 
Bubbles (2 x 10 37 erg/ s in the same energy range) , such an 
additional signal favors the AGN case, wi th th e Bubbles 
coming from the decelerated jet material [112 1. 

Yet the two signals, i.e the combined 7-ray cocoon and 
jets and the Fermi Bubbles may be created at a differ- 
ent time. Such a case would allow for a combination of 
sources, DM & AGN, Hadronic model & AGN to account 
for the total 7-ray Bubbles/haze & jets and cocoon sig- 
nals. For the hybrid scenario of DM & AGN, the edge 
at / > (left) of the south bubble could just be the 
result of the presence (overlap on the sky) of the AGN 
cocoon and thus morphologically would not need to be 
explained by DM annihilation. Similarly a cocoon in the 
north hemisphere (not claimed to be clearly revealed yet) 
could account for the north to the right edge (at Z < 0) in 
the data. Addition ally the AGN responsible for the jets 
and the cocoon of |ll2] could evacuate the cavity of the 
Bubbles which then the high energy from annihilating 
DM fill Q. 

Since the luminosity of the jets alone is only 1% of 
that of the Bubbles/haze, assuming that the current 
state of t he jets is representative of the time- averaged 
state[130], the predictions for the neutri no fl uxes from 
the DM/Hadronic cases remain the same [131 1. 
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